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We have calculated the nuclear symmetry energy Esym{p) up to densities of 4 ~ 5po with the 
effects from the Brown-Rho (BR) and Ericson scahngs for the in-medium mesons included. Using 
the Viow-k low-momentum interaction with and without such scalings, the equations of state (EOS) 
of symmetric and asymmetric nuclear matter have been calculated using a ring-diagarm formalism 
where the particle-particle-hole-hole ring diagrams are included to all orders. The EOS for symmetric 
nuclear matter and neutron matter obtained with linear BR scaling are both overly stiff compared 
with the empirical constraints of Danielewicz et al. In contrast, satisfactory results are obtained 
by either using the nonlinear Ericson scaling or by adding a Skyrme-type three-nucleon force (TNF) 
to the unsealed Viow-k interaction. Our results for Esym{p) obtained with the nonlinear Ericson 
scaling are in good agreement with the empirical values of Tsang et al. [tI and Li et al. , while 
those with TNF are slightly below these values. For densities below the nuclear saturation density 
po, the results of the above calculations are nearly equivalent to each other and all in satisfactory 
agreement with the empirical values. 

PACS numbers: 21.65.Jk, 21.65. Mn, 13.75. Cs 



I. INTRODUCTION 

The nuclear matter symmetry energy is an important 
as well as very interesting subject in nuclear and astro- 
nuclear physics. As reviewed extensively in the litera- 
ture it plays a crucial role in determining many 
important nuclear properties, such as the neutron skin 
of nuclear systems, structure of nuclei near the drip line, 
and neutron stars' masses and radii. It is especially of 
importance that constraints on the nuclear matter equa- 
tion of state (EOS) [9j and the density (p) dependence of 
the symmetry energy Egym {p) 0j Il3| up to p ~ 4po have 
been experimentally extracted from heavy-ion collisions, 
Po being the saturation density of symmetric nuclear 
matter. There have been a large number of theoretical 
derivations oi Egymi p) using , for example, the Brueckner 
Hartree-Fock (BHF) [iMll, Dirac BHF pjp-[ill, varia- 
tion al , 17] , relativistic mean field (RMF) [1^ and Skyrme 
HF [i^ many-body methods. The results of these theo- 
retical investigations have exhibited, however, large vari- 
ations for Esymip)- Depending on the interactions and 
many-body methods used, they can give either a 'hard' 
Esym{p)j in the sense that it increases monotonically with 
p up to ~ 5po, or a 'soft' one where Egymip) arises to a 
maximum value at p ~ 1.5po and then descends to zero 
at ~ 3po i, i]. It appears that the predicted behav- 
ior of Esym{p) may depend importantly on the nucleon- 
nucleon (NN) interactions and the many-body methods 
employed. 

In the present work, we shall calculate the nuclear 
symmetry energy using the low-momentum interaction 
Viow-k derived from realistic NN interactions Vnn using 
a renormalization group approach [2l| - [25i | . To our knowl- 
edge, this renormalized interaction has not yet been ap- 



plied to the study of Esym- As it is well known, most 
realistic V^j^ contain hard cores, or strong short-range 
repulsions. This feature makes these interactions not 
suitable for being directly used in nuclear many-body cal- 
culations; they need to be 'tamed' beforehand. For many 
years, this taming is enacted by way of the BHF theory 
where Vn n is converted into the Brueckner G-matrix. A 
complication of the G-matrix is its energy dependence 
(see e.g. [1^), making it rather inconvenient for calcula- 
tions. In the Viow-k approach, a different 'taming' pro- 
cedure is employed; it is performed by 'integrating out' 
the high-momentum components of Vatat beyond a deci- 
mation scale A. In this way, the resulting Viow-k is en- 
ergy independent. Furthermore Viow-k is nearly unique, 
namely the Viow-kS deduced from various realistic Vnn 
potentials (such as [27l - [30t ) are nearly identical to each 
other for decimation scale A ~ 2fm~^ [23l. |2^. 

Using this Viow-k interaction, we shall first calculate 
the equations of state (EOS) E{p, a) for asymmetric nu- 
clear matter, from which Esymip) can be obtained. Here 
E is the ground-state energy per nucleon and p is the 
total baryon density, a is the isospin asymmetry param- 
eter defined as a = (p„ — Pp)/p, where p„ and pp de- 
note, repspectively, the neutron and proton density and 
P = Pn + Pp- Our EOS will be calculated using a ring- 
diagram many-body method [3ll - [33j . As we shall discuss 
later, this method includes the particle-particle hole-hole 
(pphh) ring diagrams to all orders. In comparison, only 
the diagrams with two hole lines are included in the famil- 
iar HF, BHF and DBHF calculations. In other words, in 
these HF methods a closed Fermi sea is employed while in 
the ring-diagram framework the effects from the fluctua- 
tions of the Fermi sea are taken into account by including 
the pphh ring diagrams to all orders. 
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The nuclear symmetry energy Esym {p) is related to the 
asymmetric nuclear matter EOS by 

E{p, a) = E{p, a = 0) + E,y^{p)a^ + 0{a^). (1) 

The contributions from terms of orders higher than 
are usually negligibly small, as illustrated by our results 
in section III. With such contributions neglected, we have 



Esyrn{p)^E{p,l)-E{p,0). 



(2) 



Then the symmetry energy is just given by the energy dif- 
ference between neutron and symmetric nuclear matter. 
In calculating Esymip), the above EOS clearly play an 
important role. In our calculation, we shall require that 
the NN interaction and many-body methods employed 
should give satisfactory results for E{p,l) and E{p,0) 
of, respectively, neutron and symmetric nuclear matter. 
The use of Viow-k alone, however, has not been able to re- 
produce the empirical nuclear saturation properties, the 
predicted saturation density and binding energy per par- 
ticle being both too large compared with the empirical 
values of po — 0.16fm~^ and E ~ — 16MeV for symmetric 
nuclear matter (sil |32| . To improve the situation, it may 
be necessar y to include the effects from Brown-Rho (BR) 
scaling |34| - [3a | for the in-medium mesons, or a three- 
nucleon force (TNF) BR scaling is suitable only for 
the low density region; it suggests that the masses of light 
vector mesons in medium are reduced 'linearly' with the 
density. We consider here the EOS up to about ~ 5po 
and at such high density the linear BR scaling is clearly 
not applicable. In the present work we shall adopt the 
nonlinear Ericson scaling (4l| for the in-medium mesons 
and apply it to our Esym{p) calculations. The effects 
from the linear BR and nonlinear Ericson scalings on the 
nuclear EOS and symmetry energy will be studied. 

The organization of this paper is as follows. In sec- 
tion II we shall briefly describe our derivation of the low- 
momentum interaction Viow-k using a T-matrix equiv- 
alence approach. Some details about the calculation of 
the EOS for asymmetric nuclear matter from this inter- 
action with the pphh ring diagrams summed to all orders 
will also be presented. The Ericson scaling is a nonlinear 
extension of linear BR scaling. The difference between 
them will be addressed in this section. Our results will be 
presented and discussed in section III. A summary and 
concision is contained in section IV. 



II. FORMALISM 

We shall calculate Esvmip ) using a low-momentum 
ring-diagram approach [3ll433[, where the pphh ring dia- 
grams are summed to all orders within a model space of 
decimation scale A. In this approach, we employ the low- 
momentum interaction Viow-k |2ll425| . Briefly speaking, 
this interaction is obtained by solving the following T- 



matrix equivalence equations: 
T{k',k,k^) = X4^N(fc',fc) 
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VMk\q)T{q,k,k^) 
k^ — (p- 



q^dq, (3) 
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Viow-k(fc',<7)Tlow-k(g, k,P 



-q^dq, (4) 



T(fc', fc, = Tio^_k{k', k, (k', k) < A. (5) 

In the above Vnn represents a realistic NN interaction 
such as the CDBonn potential [ItJ . P denotes principal- 
value integration and the intermediate state momentum q 
is integrated from to oo for the whole-space T and from 
to A for Tiow-k- The above Viow-k preserves the low- 
energy phase shifts (up to energy A^) and the deuteron 
binding energy of Vnn- Since Viow-k is obtained by inte- 
grating out the high- momentum components of Vnn, it 
is a smooth 'tamed' potential which is suitable for being 
used directly in many-body calculations. 




FIG. 1: Sample ring diagram included in the equation of state 
E{p,a). Each wave line represents a Viom-k vertex. The HF 
one-bubble insertions to the Fermion lines are included to all 
orders. 



We use a ring-diagram method [31l - l33| to calculate the 
nuclear matter EOS. In this method, the ground-state 
energy is expressed as E{p, a) = Ef^'^'^{p, a) -f AE{p, a) 
where E^'^^'^ denotes the free (non-interacting) EOS and 
AE is the energy shift due to the NN interaction. In 
our ring-diagram approach, AE is is given by the all- 
order sum of the pphh ring diagrams as illustrated in Fig. 
1. Note that we include three types of ring diagrams, 
the proton-proton, neutron-neutron and proton-neutron 
ones. The proton and neutron Fermi momenta are, re- 
spectively, kpp = (37r^pp)^/*^ and kpn = (37r^p„)^/^. 
With such ring diagrams summed to all orders, we have 
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m ijkl<A 

xy;;(fc?,A)(ij|Mow-k|A:/), (6) 

where the transition amplitudes Y are obtaind from a 
pp/i/i RPA equation [3ll - l33l |. Note that A is a strength 
parameter, integrated from to 1. The above ring- 
diagram method reduces to the usual HF method if 
only the first-order ring diagram is included. In this 
case, the above energy shift becomes AE{p,a)HF = 
^Y^riinj{ij\Viow-]<i\ij) where nfe=(l,0) if fc(<, >)kFp for 
proton and nfc=(l,0) if fc(<, >)kFn for neutron. 

It is well known that the use of the free-space Vnn 
alone is not adequate for describing nuclear properties at 
high densities. To satisfactorily describe such properties, 
one may need to include the three-nucleon force [40 ] or 
the in-medium modifications to the nuclear interaction. 
In the present work, we shall employ in our EOS calcu- 
lations nuclear interactions which contain the in-medium 
modifications suggested by the Brown- Rho (BR) [sj, HH] 
and Ericson l4ll sca ling s. These scalings are based on 
the relation j34 Issl . l39l | that hadron masses scale with 
the quark condensate (qq) in medium as 
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where m* is the hadron mass in a medium of density p, 
and m is that in free space. The quark condenstate {qq) 
measures the chiral symmetry breaking, and its density 
dependence in the low-density limit is related [13, IS^I to 
the free ttN sigma term SjrAr by 



{mm 



= 1 - 



N 



(8) 



where /tt = 93MeV is the pion decay constant and SttTV = 
45 ± 7MeV [l^j. Applying the above scaling to mesons 
in low-density nuclear medium, one has the linear scaling 

M 



m po 



(9) 



where m* and m are, respectively, the in-medium and 
free meson mass, and C is a constant of value ^ 0.15. 
The above scaling will be referred to as the linear BR scal- 
ing. Nucleon-nucleon interactions are mediated by meson 
exchanges, and clearly the in-medium modifications of 
meson masses can significantly alter the NN interaction. 
These modifications could arise from the partial restora- 
tion of chiral symmetry at finite density/temperature 
or from traditional many-body effects. Particularly im- 
portant are the vector mesons, for which there is now 
evidence from both theory [35l [isl l46j and experiment 
[47I l48j that the masses may decrease by approximately 



10 — 15% at normal nuclear matter density and zero tem- 
perature. (Pions are not scaled because they are pro- 
tected by chiral symmetry.) Density-dependent nuclear 
interactions obtained by applying the above scaling to 
the light mesons (w, p and a) which mediate the NN 
potential have been employed in studying the properties 
of nuclear matter [HI, IM IH, [13 and the ^^C N 
P-decay 

We are interested in the EOS and Egym up to densities 
as high as p ~ 5po, and at such high densities the above 
linear scaling is clearly not suitable. How to scale the 
mesons in such high density region is still by and large 
uncertain. We shall adopt here the Ericson scaling 41 1 
which is an extension of the BR scaling. In this scaling, 
a new relation for the quark condensate (qq) based on 
chiral symmetry breaking is employed, namely 
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Note that this relation agrees with the linear scaling re- 
lation of Eq.(8) for small p. The above scaling suggests 
a non-linear scaling for meson mass 
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(11) 



with D — ^jTj^i , and we shall refer to this scaling as 
the nonlinear Brown-Rho-Ericson (BRE) scaling. Using 
the empirical values for {Ti-kN, Po, /tt, "t-tt), we have D= 
0.35±0.06. In the present work, we shall employ the one- 
boson exchange BonnA potential with its (p, cj, a) 
mesons scaled using both the linear (Eq.(9)) and nonlin- 
ear (Eq.(ll)) scalings. This potential is chosen because 
it has a relatively simple structure which is convenient 
for scaling its meson parameters. 



III. RESULTS AND DISCUSSIONS 

Using both the unsealed and scaled BonnA potentials, 
we first calculate the ring-diagram EOS for symmetric 
nuclear matter to investigate if they can give saturation 
properties in good agreement with the empirical values. 
We employ the low-momentum interactions Viow-k from 
these potentials using a decimation A — 3.0/to~^, which 
is chosen because we are to study the EOS up to high 
densities of ^ 5po. As shown in Fig. 2, the EOS (la- 
belled 'V/oui-fe alone') calculated with the unsealed po- 
tential saturates at kp ~ 1.8/m^^, which is too large 
compared with the empirical value, and it also overbinds 
nuclear matter. We then repeat the calculation includ- 
ing the medium modifications from the BR scalings. For 
the linear BR scaling (Eq.(9)), we have used Ctj=0.128, 
Cp=0.113 and Cct=0.102. These parameters are chosen 
so as to have satisfactory saturation properties, namely 
they give Eq/A ~-15.5 MeV and po ~ 0.17 /to'^. in Fig. 
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2 we also present our results obtained with the nonlinear 
BRE scahng (Eq.(ll)) using parameters = Dp=Q.AQ 
and -Dcr=0.30. They were chosen to provide satisfactory 
results for Eq/A and po- It is of interest that for den- 
sities < po the EOS given by the linear BR and nonlin- 
ear BRE scalings are practically eqivalent to each other. 
From Eqs.(7-ll), we see that the parameters D and C 
obtained from the density dependence of quark conden- 
sates is 0.29 <D< 0.41 and C ~ D/i. It is noteworthy 
that the C and D parameters we have employed in the 
EOS calculations agree well with the above theoretical 
values. 

As also seen from Fig. 2, the above equivalence be- 
gins to disappear for densities larger than po- There the 
EOS given by the linear scaling is much stifFer than that 
given by the nonlinear one; the difference between them 
becomes larger and larger as density increases. In ad- 
dition to the above two EOS, we have also calculated 
an EOS using the interaction given by the sum of the 
unsealed Viow-k and an empiriral Skyrme three-nucleon 
force (TNF). The well-known emipirical Skyrme force 
I49II is of the form 



i<j i<j<k 



(12) 



where V{i,j) is a two-nucleon momentum dependent in- 
teraction, and V(i,j,k) is a zero-range three-nucleon inter- 
action which has played an indispensible role for nuclear 
saturation. For nucleons in a nuclear medium of density 
p, this three-nucleon force becomes a density-dependent 
two-nucleon force commonly written as 



Vp{i,j) = —pSin 



rl)- 



(13) 



In Fig. 2 the EOS labelled 'V/o™_fc with TNF' is obtained 
using the combined interaction of Viow-k (unsealed) and 
Vp. The parameter is adjusted so that the resulting 
EOS gives satisfactory saturation properties for symmet- 
ric nuclear matter. The EOS shown has <3=2000 MeV- 

It is of interest that the above three EOS (linear BR 
and nonlinear BRE, TNF) are nearly identical for densi- 
ties < po, but they deviate from each other with increas- 
ing densities. Without experimental guidelines about the 
nuclear matter EOS above po, it would be difficult to de- 
termine which of these three EOS has the correct high 
density behavior. Fortunately, heavy-ion collision exper- 
iments conducted during the last several years have pro- 
vided us with constraints of the EOS at high densities. 
Danielewicz et al. [9J have obtained a constraint on the 
EOS for symmetric nuclear matter of densities between 
2po and 4.5poj as shown by the red solid-line box in Fig. 
3. Comparing our three EOSs with their constraint, the 
linear BR EOS is clearly not consistent with the con- 
straint and should be ruled out. This linear scaling is 
suitable for low densities, but definitely needs modifica- 
tion at high densities. It is primarily for this purpose that 
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FIG. 2: Ring-diagram EOSs calculated with the Viow-k inter- 
action alone, with the linear BR scaling of Eq.(9), with the 
nonlinear BRE scaling of Eq.(ll), and with the addition of a 
Skyrme-type three-nucleon force (TNF) . 



we have considered the nonlinear scaling. As displayed 
in Fig. 3, the EOS with the nonlinear BRE scaling is in 
much better agreement with the constraint than the lin- 
ear BR one. It satisfies the constraint well except being 
slightly above the constraint at densiies near ~ 4.5po- It 
is of interest that the EOS using Viow-k with the Skyrme- 
type TNF exhibits even better agreement with the con- 
straint. 
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FIG. 3: Comparison of the calculated equations of state for 
symmetric nuclear matter with the constraint (solid-line box) 
of Danielewicz See text for more explanations. 

So far we have studied the effects of the BR scalings 
and the TNF three-nucleon force on the EOS for symmet- 
ric nuclear matter. The neutron matter EOS is also an 
interesting and important topic ^50, 51]. It plays a crucial 
role in determining the nuclear symmetry energies as well 
as the properties of neutron stars. It should be of interest 
to study also the effects of the above BR/BRE scalings 
and the TNF force on the EOS of neutron matter. Us- 
ing the same Viow-k ring-diagram framework employed 
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for symmetric nuclear matter and the same C, D and 
parameters, we have caculated the neutron matter EOS 
up to 4.5po- Our calculated neutron-matter EOS are dis- 
played in Fig. 4. Danielewicz et al. ^ have given two 
different constraints for the neutron matter EOS: a stiff 
one (upper black solid-line box) and a soft one (lower red 
solid-line box) which are both displayed in Fig. 4. As we 
can see, the linear BR EOS is again producing too much 
pressure. The nonlinear BRE EOS agrees well with the 
stiff constraint (upper box) while the TNF EOS is fully 
within the soft constraint box. To further test these two 
EOS (nonlinear BRE and TNF), it would be very helpful 
to have narrower experimental constraints on the neutron 
matter EOS. 
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FIG. 4: Comparison of the calculated equations of state for 
neutron- matter with the constraints of Danielewicz See 
text for more explanations. 

The symmetry energy Esym is a topic of much current 
interest, and extensive studies have been carried out to 
extract its density dependence from heavy-ion collision 
experiments 0, E^- Based on such experiments, Li et al. 
suggested an empirical relation 

Esym{p) ~ 31.6(p/po)^; 7 - 0.69 - 1.1, (14) 



exp' 

m 



for constraining the density dependence of the symme- 
try energy. Also based on such experiments, Tsang et 
al. 0] recently proposed a new empirical relation for the 
symmetry energy, namely 



Esym (p) 



s.k 



2/3 



C. 



(15) 



where Cs,k = 25MeV, Cs^p = 35.2MeV and 7, « 0.7. It 
should be useful and of interest to check if our calculated 
Esymip) is consistent with the above relations. 

Using the ring-diagram framework described earlier, 
we have calculated the ground-state energy E{p, a) for 
asymmetric nuclear matter. (Recall that the asymmetry 
parameter is a = (/9„ — Pp)/p.) Some representative re- 
sults are shown in Fig. 5: the results in the left panel are 



FIG. 5: Ring-diagarm equations of state of asymmetric nu- 
clear matter. See text for more explanations. 



obtained with the ^Viow-k with TNF' interaction while 
for the right panel the 'nonlinear BRE' interaction is 
used. As seen, E(p, a) varies with almost perfectly 
linearly, for a wide range of p. (Note that in Fig. 5 
we plot the energy difference Esym{p,o.) — Esymip,0)-) 
This is a desirable and remarkable result, indicating that 
our ring-diagram symmetry energy can be accurately ob- 
tained from the simple relation given by Eq.(2), namely 
the energy difference between neutron and symmetric nu- 
clear matter. 

In Fig. 6, the 'shaded area' represents the empirical 
constraint, Eq.(14), of Li et al. [l3- As seen, there are 
large uncertainties in the high-density region. The em- 
pirical relation Eq.(15) of Tsang et al. [7| is given by 
the 'second curve from bottom' in the figure. As seen, 
the density dependence of this relation is slightly below 
the softest limit (lower boundary of the shaded area) of 
Eq.(14). Our 'nonlinear BRE' results are in the middle 
of the shaded area, in good agreement with the empirical 
constraint of [l^l- Our results with the TNF force are 
below the empirical ones of both 

and 0, giving a 
softer density dependence than both. It may be noticed 
that for densities (p < po), the calculated and empirical 
results are all in good agreement with each other. The 
symmetry energies given by them at po are all close to 
~ 30MeV, which is also the only well determined em- 
pirical value. Furthermore, our calculated symmetry en- 
ergies all increase monotonically with density. We have 
required our nuclear matter EOS to satisfy certain em- 
pirical constraints, and with such requirements it may 
be difficult for our present calculations to have a soft 
Esym{p) as soft as the supersoft one of Q which satu- 
rates at density near ^ 1.5po- 

We have found that our symmetry energies can be well 
fitted by expressions of the same forms as Eqs.(14) and 
(15), with the exponents 7 and 7^ treated as parameters. 
In Table I, we compare the exponents determined from 
our results with the empirical ones of jl3l and Q . The 7 
exponent given by the nonlinear BRE scaling is in good 
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agreement with the empirical values of [lo!]. The empir- 
ical 7i of is, however, about half-way between the 7^ 
obtained with 'nonlinear BRE' and that with 'TNF'. 



200 




FIG. 6: Comparison of the density dependence of our cal- 
culated symmetry energies with the empirical results of 7] 
(dot-dash line) and tlQi] (shaded area). 



TABLE I: Comparison of the density exponents for the nu- 
clear symmetry energy Esym{p)- The exponents 7 and 7i are 
defined respectively in Eqs. (14) and (15). 
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IV. SUMMARY AND CONCLUSION 

Employing the Viow~k low-momentum interactions, we 
have calculated the nuclear symmetry energy Esym (p) up 
to a density of ~ 5po using a ring-diagram framework 
where pphh ring diagrams are summed to all orders. We 
first calculate the EOS for symmetric nuclear matter and 
neutron matter and compare our results with the corre- 
sponding empirical constraints of Danielewicz et al. Q . 
To have satisfactory agreements with such constraints, 
we have found it necessary to include certain medium cor- 
rections to the free-space NN interations. In other words, 
the effective NN interactions in medium are different 
from those in free space, and when using them in nuclear 
many-body problems it may be necessary to include the 
renormalization effects due to the presence of other nucle- 
ons. We have considered several methods to incorporate 
such medium corrections. Although the nuclear matter 



saturation properties can satisfactorily be reproduced by 
including the medium corrections from the well-known 
linear Brown-Rho scaling for the in-medium mesons, this 
scaling produces an EOS which is too stiff compared 
with the Danielewicz constraints. We have found that 
the EOS obtained with the nonlinear Brown- Rho-Ericson 
scaling are in good agreement with the Danielewicz con- 
straints. We have considered another method to ren- 
der the effective interaction density dependent, namely 
adding a Skyrme-type three-nucleon force (TNF) to the 
unsealed Viow-k interaction. The EOS so obtained are 
also in good agreement with the Danielewicz constraints, 
but the resulting neutron matter EOS is significantly 
softer than that with the nonlinear scaling. The three 
methods (linear and nonlinear scalings, and TNF) all 
have reproduced well the empirical saturation properties 
of nuclear matter {po « O.lTfm"^ a.ndEo/A w -15MeV), 
but their results at high densities are different. We have 
determined the scaling parameters C (linear BR scaling) 
and D (nonlinear BRE scaling) by fitting the above sat- 
uration properties. It is encouraging that the results, 
(0.102 < C < 0.128) and (0.30 <D< 0.40), so obtained 
are actually in good agreement with the theoretical result 
D ~ 0.35 ± 0.06 ~ 3C given by Eqs.(7-ll). 

Including the above medium modifications, we pro- 
ceed to calculate the nuclear symmetry energies. We 
have found that the Esym{p, a) given by our asymmet- 
ric ring-diagram calculations depends on almost per- 
fectly linearly. This is a rather surprising and useful re- 
sult, suggesting that the symmetry energy can be reli- 
ably obtained from the simple energy difference between 
symmetric nuclear matter and neutron matter. Our sym- 
metry energies obtained with the nonlinear BRE scaling 
agree well with the empirical constraints of (lo| . and are 
slightly above the empirical values of Q. Our results 
with the TNF force is slightly below the empirical re- 
sults of both [l3l and Q . The non-linear Ericson scaling 
has given satisfactory results for the equations of states 
of nuclear matter and nuclear symmetry energies up to 
a density of ^ 5po- We believe this scaling provides a 
suitable extension of the linear BR scaling to moderately 
high densities of < 5po- Our calculated Esym{p) a-U in- 
crease monotonically with p up to ^ 5po- It may be of 
interest to carry out further studies about the possibil- 
ity of obtaining a supersoft symmetry energy which may 
saturate at some low density of ~ l.Spo @- 
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